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Abstract 

A strongly correlated electron system associated with the quantum superalgebra 
U 9 [osp(2|2)] is studied in the framework of the quantum inverse scattering method. 
By solving the graded reflection equation, two classes of boundary-reflection K- 
matrices leading to four kinds of possible boundary interaction terms are found. 
Performing the algebraic Bethe ansatz, we diagonalize the two-level transfer ma- 
trices which characterize the charge and the spin degrees of freedom, respectively. 
The Bethe-ansatz equations, the eigenvalues of the transfer matrices and the energy 
spectrum are presented explicitly. We also construct two impurities coupled to the 
boundaries. In the thermodynamic limit, the ground state properties and impurity 
effects are discussed. 
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1 Introduction 



In recent decades, there has been considerable interest in the study of strongly 
correlated electron systems in reduced dimensions exhibiting non-Fermi-liquid 
behaviour [1-3] . This has been motivated by the surprising properties of high- 
T c superconductors [4], heavy- fermion alloys and compounds [5], and the con- 
cept of Luttinger liquids [6-8] . Similar non-standard behaviour that lies out- 
side the realm of Fermi liquid theory was also observed in the magnetic prop- 
erties of systems displaying the Kondo effect [9]. The Kondo problem was 
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solved exactly by means of the Bethe ansatz and the quantum inverse scat- 
tering method (QISM) [10]. Following this approach, many one-dimensional 
integrable quantum systems with impurities [11] have been constructed as 
inhomogeneous solutions of the Yang-Baxter equations (YBE) [12]. 

A particularly intriguing situation corresponds to perfect backscattering impu- 
rities which can be realized conveniently by integrable open boundary condi- 
tions (BC). A systematic approach to handle quantum systems with backscat- 
tering boundaries is provided by Sklyanin's work on the reflection equations 
(RE) [13]. In analogy with the YBE in the bulk, the RE guarantee the fac- 
torization of the scattering matrices at the boundaries. Thus, starting from a 
solution of the YBE, which yields a solvable model with periodic BC, one con- 
structs suitable integrable boundary conditions such that the RE are fulfilled 
[14-17]. Due to integrability, the bulk impurities obtained by inhomogeneous 
solutions of the YBE are pure forward scatterers. Thus their combination with 
the backscattering boundaries may be expected to model physically relevant 
impurity systems. In the context of boundary integrable quantum field theories 
[18], a model with impurities can be mapped onto a model with certain bound- 
ary conditions. Because impurity effects play a decisive role for the transport 
in quantum wires, general boundary conditions for strongly correlated elec- 
tron systems open many opportunities to investigate the transport properties 
in Luttinger liquids or quantum wires. A proper boundary field may have a 
feasible realization by applying boundary external voltages and external mag- 
netic fields in experiments on quantum wires [19]. Moreover, we may expect 
that the local state induced by the boundary fields inherits signatures of the 
bulk Luttinger liquid [6-8] . The physical quantities such as magnetization, the 
compressibility, susceptibility, and specific heat, etc., may be manipulated by 
the Bethe-ansatz equations. For this reason, integrable models that combine 
bulk and boundary impurities have recently attracted much attention [20-23] . 

Of particular interest are strongly correlated electron systems associated with 
super symmetric solutions of the YBE. Models corresponding to non-excep- 
tional Lie superalgebras as, for instance, gl(2|l) and osp(2|2) [24-26], have 
provided interesting non-perturbative information [27,28] for generalizations 
of well-known models such as, e.g., the Hubbard model [29]. A further gener- 
alization was achieved by considering the solution of the YBE related to the 
quantum superalgebra U g [osp(2|2)] [25]. This model has two fermionic and two 
bosonic degrees of freedom. Furthermore, the Hamiltonian [26] corresponds to 
a lattice regularization of the integrable double sine-Gordon model [30] . The 
continuum version of this model with boundary fields is known to describe 
tunneling effects in quantum wires [31]. These wires are believed to represent 
a realization of Luttinger liquids. We remark that the coordinate Bethe ansatz 
for an open U g [osp(2|2)] chain has recently been studied [32]. 

In the present work, we perform the algebraic Bethe ansatz for a supersymmet- 
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ric open chain associated with the quantum superalgebra U 9 [osp(2|2)]. Due to 
the super symmetric structure of the model, we use the graded version of the 
QISM [28]. Starting from the U g [osp(2|2)] solution of the graded YBE, we con- 
struct a supersymmetric correlated electron system with boundary fields by 
solving the graded RE. We find that the model contains a hidden anisotropic 
XXZ open chain characterizing the spin degrees of freedom [1]. This plays 
a crucial role in our solution, which proceeds in two steps, treating charge 
and spin degrees of freedom separately by a nested graded Bethe ansatz. This 
structure also suggests a natural way to incorporate impurities coupling to 
charge and spin degrees of freedom while preserving the integrability of the 
model. From the Bethe-ansatz equations, we obtain the ground state energy at 
half-filling in the thermodynamic limit. We also discuss integrable impurities 
coupled to the boundaries. 



The paper is organized as follows. In section 2, we present the U 9 [osp(2|2)] 
solution of the YBE and solve the corresponding RE. Furthermore, we give 
an explicit expression of the Hamiltonian in terms of fermionic operators. 
Section 3 is devoted to the derivation and the solution of the Bethe-ansatz 
equations by means of the graded QISM. In section 4, the ground state energy 
in the thermodynamic limit is obtained for the model with boundary fields 
and impurities. We conclude in section 5. 



2 The solutions of the graded RE 



We begin by considering a two dimensional classical lattice model, where, to 
each bond of the lattice, two bosonic and two fermionic degrees of freedom are 
associated. It is a vertex model, hence interaction takes place at each vertex, 
and the energies of the various local configurations determine the statistical 
weight of the configuration. The corresponding Boltzmann weights form a 
2 4 x 2 4 quantum i?-matrix with grading 'bffb', where 'b' stands for bosonic 
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and 'f ' for fermionic. It has the form 
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x — q 1 — q q(x — 1) 

Wi = - J, W 2 = - J, W 3 = 

1 — xq z 1 — xq A 1 — xq 
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2w 2 1 — x 

W 4 = XW 2 , W 5 = — , Wq = — W 2 , 

1 + X 1 + X 

1 + q 2 x 2xw 2 x + q 2 

w 7 = — r w 6 , w 8 = -— — , w g = — -Wa, 

1 — q A 1 + x 1 — q z 

2x 2 w 2 

w w = xw 6 , w u = — , (2) 

1 + x 

with x = e A characterizing the difference A of the pseudo-momenta of the 
particles whose two-body scattering is described by the quantum i?-matrix 
(1); q is the deformation parameter. 

In what follows, we are going to use also the standard diagrammatic repre- 
sentation [10], where the i?-matrix corresponds to the two-particle scattering 
picture 



R12W = X • ( 3 ) 




The quantum i?-matrix (1) satisfies the graded YBE 

i? 12 (A - z/) J R 13 (A) J R 23 (z/) = #23(^i3(A)i?i 2 (A - v) (4) 



guaranteeing the integrability of the model with periodic BC, i.e., the fac- 
torization of the scattering matrices into two-body scattering matrices [12]. 
Above, Rij with % — 1,2 and j = 2,3, denotes on which of the ith and jth 
spaces of V\ ® s V 2 ® s V 3 the i?-matrix acts. In the remaining space, acts 
an identity. Here ®s denotes the graded tensor product 

[A® s BU, l5 = {-l)^ +p ^ p ^A ai B p5 , (5) 

with the Grassmann parities obey P(l) = -P(4) = and P(2) = P(3) = 1 
with respect to the grading 'bffb'. 

For other BC such as twisted and open BC, the graded YBE will still account 
for the bulk part of the model, but the boundary terms have to be chosen 
appropriately in order to preserve the integrability of the model. In particu- 
lar, solutions of the RE yield integrable vertex models with open (reflecting) 
boundaries or the equivalent integrable quantum spin chains with boundary 
fields. In the context of two-body scattering, the RE characterize the consis- 
tency conditions for the factorizations of the two-body boundary scattering 
matrices at boundaries. Taking into account the grading 'bffb', the left and 
right reflection matrices, K_ and K + , are required to satisfy the following RE 

Rn(\ - v) K- (A) R 21 (X + u) K- {y) 

= K-{v)Ri*{\ + v) k-(\)R 2 i(\-v), (6) 

Rf^(v - A) ir^(A) Rfl^i-x - u) k^ 2 (v) 

= kf 2 (u)RH 1 ^ t2 (-\-u) k + \X) Rfl^iv - A), (7) 
respectively. Here, we used the conventional notation 

X= X ® s Iy 2 , X= I Vl ®5 X, (8) 

where Iy denotes the identity operator on V, and, as usual, R 2 i = P ■ Ru • P- 
Here P is the graded permutation operator which can be represented by a 
2 4 x 2 4 matrix, i.e., 

P a ^s = (-l) P{a)PW S aS 6^. (9) 

Furthermore, superscripts St a and St a denote the supertransposition in the 
space with index a and its inverse, respectively, 

(^■) St = (-ljlWl^, (Aijf = (-lyW+WWAji. (10) 
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We found two solutions of the RE (6) and (7) for diagonal boundary K±- 
matrices 



(K W (X) 
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'ID 



K + (X) 



K^(X — > —A — i7r, Z-^-i-). 



(12) 



In the first solution, the entries are given as 



K ^(X) = (x + ^q)(x-^q- 1 ), 

K^(X) = K^(X) = (x-i+t-qXx-t-q- 1 ), 

K ^(X) = (x- 1 +^q)(x~ 1 -^q- 1 ), (13) 
whereas the second solutions corresponds to 



(x + £~q)(x - C 1 ?)^ - C-'T 1 ), 
(x^ + Z-q^x-C^Kx-Z-q- 1 ), 

(x + ^(x- 1 - C 1 ?)^ - Z-q- 1 ), 

(X + C q)(x ; - C'q)^- 1 - i q )■ (14) 

Here x = e A as defined in (2) and £± are free parameters characterizing 
the boundary fields and the boundary interactions. We emphasize that the 
permutation-transposition symmetry 

i?st 1 st 2(A) = PjRi2(A)P) (15) 



K W (X) = 
^i 2) (A) = 
^i 3) (A) = 
K^(X) = 



the unitarity 

R 12 (X)R 21 (-X)=1, (16) 
and the graded crossing symmetry 



i2 12 (A) = VR S £(-X - in) V- 1 Ci c A (A) i?r) , (17) 
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with 
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[l-xq 2 ) 
(l-x) ' 



(18) 



not only result in the isomorphism (12) between K + (X) and K_(\), but also 
constitute the necessary ingredients for the integrability of the model with 
boundaries. We also remark that the grading does not play a role for the di- 
agonal K_(A)-matrix (13), which coincides with the one for the non-graded 
i?-matrix given in Ref. [32]. However, its companion K + (X) given by (12) is 
different from that in the non-graded case because it has to obey the graded 
crossing symmetry (17). It is worth emphasizing that the boundary parame- 
ters £_ and should inherit the same crossing property as that imposed on 
the pseudo-momenta A such that the boundary terms of the corresponding 
Hamiltonian are completely symmetric. In Appendix A we present an ansatz 
to work out the solutions (13) as well as (14). 

A more important object in the context of the QISM is the transfer matrix 
of an integrable system, which can be considered as a generating function of 
the infinite integrals of motion due to its commutativity for different values of 
the spectral parameter. Actually, the RE (6) and (7) together with the YBE 
(4) and the symmetries of the quantum i?-matrix ensure the commutativity 
of the double-row transfer matrix 



1 2 



L 



r(A) 



Str \k + (X)T(X)K4X)T- 1 (-X) 



(19) 



Here Stro denotes the supertrace carried out in auxiliary space Vq. The mon- 
odromy matrices T(A) and T _1 (— A) are defined by 



T(A) 



1 2 



— Rl,o(X)Rl-1,o(X) ■ ■ ■ R2fl(^)Rlfl(ty, 



(20) 
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12 L 



T-\-\) = 



(21) 



respectively. The Hamiltonian associated with the quantum i?-matrix (1) is 
related to the double- row transfer matrix (19) as 



r(A) = ciA + c 2 {H + const. )A 2 + . . . , 



(22) 



where Cj, % — 1, 2, are scalar functions of the boundary parameters £±. Taking 
into account the Grassmann parity of the host vertices, after some lengthy 
algebra, one can present the Hamiltonian explicitly in terms of the fermionic 
creation and annihilation operators c] CT and Cj a acting on the site j and car- 
rying the spin index a — ±. It is given by 



L-l 



H = E H iJ+i + U(S-)n 1>+ n 1 ,- + U(S+)n L>+ n L ,. 
j'=i 



(23) 



cr=± 



where n Ji0 - = c] j(7 c ■ is the fermion number operator, and the bulk Hamiltonian 
is chosen as 



<T=± 

-(<? - S [ c L c i+i,<r - hx -] - n i>-^] 

+(? + [4+4-^+1 - c i+i.+ ~ c ],+ c l+i ,- c i+i,+ c i - + h - c - 



+ (<? + <? *) + _ + n ji+ n j+h 

+2(q + q- 1 % 



(24) 



with rij = nj t+ + rij t - and I the identity operator. The remaining terms in (23) 
are the on-site Coulomb coupling t/(£±) and the chemical potentials Vi j(r and 
Vl j(T at the ends of the chain. For the first solution (13), we find 



2(q 2 - q- 2 )qi ± 

(?-e±)(i+?e±)' 

(g-jT%f + e-) 



(9 - *+) ' 



(25) 
(26) 
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whereas the second solution (14) yields the same expression (25) for U(£±), 
Vi _ and Vt _ but 

v fa - r'M- - 1) „ fa -<r')fa; + -i) ,„ 7 , 



as presented in Appendices A and B. The Hamiltonian (23) contains hop- 
ping terms with occupation numbers, double hopping terms, on- and off-site 
Coulomb interaction in the bulk, as well as on-site Coulomb interaction and 
chemical potentials at the boundaries. We notice that the boundary terms 
corresponding to the first solution (13) act as boundary chemical potentials, 
whereas in the second case (14) they act as boundary magnetic fields. The 
two cases could provide four possible classes of boundary conditions, which 
lead to different boundary shift factors in the Bethe-ansatz equations of the 
model. In order to keep the Hamiltonian hermitian, we restrict ourselves to 
q = e 17 and boundary parameters £± = e^ with real 7 and £ ± . We note that 
the Hamiltonian in Ref. [26] differs from the bulk part (24), but that they are 
related to each other by a canonical transformation. Nevertheless, it will be 
shown that this transformation does not change the Bethe-ansatz equations 
in the bulk. 

So far, we finished the first step towards the algebraic Bethe-ansatz solution for 
the model associated with the quantum i?-matrix (1). Next we shall proceed 
with the factorization of the transfer matrix (19). In this paper, we restrict 
ourselves to the first solution (13) of the RE, which leads to a perfect factor- 
ization of the transfer matrix (19) acting on the pseudo-vacuum state. The 
second class of solutions of the RE, given in (14), leads to a very complicated 
factorization of the transfer matrix (19) acting on the pseudo- vacuum state 
and the multi-particle states. However the ansatz formulated in the following 
section works in a similar way for the model with other boundary terms. 



3 The algebraic Bethe-ansatz approach 



In order to accomplish the algebraic Bethe ansatz for an integrable system 
with boundaries, we first need to diagonalize the transfer matrix of the model 
acting on the pseudo-vacuum state. As usual, we rewrite the transfer matrix 
(19) in the following form 

r(A) = Str [^+(A)C/_(A)], (28) 
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where U-(X) is defined by 



1 2 



C/_(A) = 



T(A)K_(A)T- 1 (-A). (29) 



One can verify that U-(X) also satisfies the graded RE (6). With regard to the 
structure of the i?-matrix (1), we choose the standard ferromagnetic pseudo- 
vacuum state [34] 

|0) = |0) L <8> • • • <8> |0)i (8) • • • (8) |0)i, (30) 



where |0)j = (^^j ® acts as a highest-weight vector. Following Refs. 

[34,35], we label the elements of the monodromy matrix T(A) by 



T(A) = 



( B{X) B 1 (X) B 2 (X) F(X)\ 

d(A) A n (X) A 12 (X) E 1 (X) 

C 2 (A) A 21 (A) A 22 (A) £ 2 (A) 

VC 3 (A) C 4 (A) C 5 (A) -D(A) / 



(31) 



and further 



-A) = 



C/_(A) 



/S(A) 


B 1 (X) 


B 2 (X) 


F(X)\ 


C 1 (X) 


A~n(X) 


A 12 (X) 


^i(A) 


C 2 (X) 
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\Cs(X) 


64(A) 


C 5 (X) 


D(X)J 


(B{X) 
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B 2 (X) 


m\ 


C 1 (X) 


A 11 (X) 


A 12 (X) 


^i(A) 


C 2 (X) 


A 21 (X) 


A 22 (X) 


E 2 (X) 


V 63(A) 


C 4 (X) 


65(A) 


D(X)J 



(32) 



(33) 



From the structure of the i?-matrix (1), the relation (16) — which one uses to 
construct the inverse i?-matrix — and the definitions (20) and (21), one can 
deduce that the operators B a (X) and B a (X) (a = 1,2) act as creation fields 
acting on the reference state, creating particles with pseudo-momenta A and 
—A, respectively. While E a (X) and E a (X) are the 'dual' creation fields to -B a (A) 
and B a (X), the operators Cj(A) and Ci(X) (i = 1, • • • , 5) behave as annihilation 
fields. Furthermore, using an invariant of the Yang-Baxter algebra 
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-A)i? 12 (2A) T(A) =T(A)i? 12 (2A) T^-A), 



(34) 
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we obtain, apart from an overall factor Q(X) = K ( l\X)K ( ^\X), the eigenvalue 
of the transfer matrix as 

r(A)|0> = {w?(A)5(A) + £ W^ +1 (A)i aa (A) + W+(A)£(A)j |0), (35) 
where we introduced the transformations 



A aa (X) = A aa (X) - ^—^B(X) = W a+1 (X)A aa (X)A aa (X), (36) 
= W^(X)D(X)D(X). (37) 



Here, 



Wr(X) = l, (38) 

W +(\) - ^ + ~ 1 ^ X + - £+) r . q x 

1 1 J (^ 2 - <? 2 x* 2 + i)(x - tf+xs? + e + ) ' 1 ' 

w + (\)- x2 ( x <iZ+- 1 )( x Z+ + q) (A2 \ 

= (a: 2 ~ 1)QeV + ~ 1)«- + g) ^ 

41 ' x 2 (x 2 + 1)(iV -l)(a;g-f_)(a; + g£_)' 1 J 

5(A)|0)= £l (A)|0) = H/r(AX(A)|0), (44) 

iaa(A)|0) =e a+ i(A)|0) = W- +1 (X)w 2 3 L (X)\0), (45) 

D(A)|0) =e 4 (A)|0) = «; 4 -(A)«; 2L (A)|0). (46) 

In Appendix C, we give some useful relations for the eigenvalue problem (35), 
the factorization of the transfer matrix (28) acting on the pseudo- vacuum state 
(30). We note that the operators B a (X), a = 1,2, constitute a two-component 
vector with both positive and negative pseudo-momenta still playing the role 
of the creation fields acting on the pseudo- vacuum state. In Eq. (33), E a (X) 
are the components of the dual creation fields, whereas one can show that the 
operators Cj(A) are still the annihilation fields acting on the pseudo- vacuum 
state. The integrability of the model leads to a perfect factorization of the 
transfer matrix acting not only on the pseudo-vacuum state but also on the 
multi-particle states. Conversely, the factorization of the transfer matrix on the 



11 



pseudo-vacuum state reveals the consistency between the boundary reflection 
-fT±-matrices and the integrability of the model. 

In order to make further progress we return to the graded RE (6) and derive 
commutation relations between the diagonal fields and the creation fields. 
In general, the algebraic Bethe-ansatz solution for an integrable model with 
open BC is more complicated than that in the case of periodic BC due to 
the appearance of positive and negative rapidities. As in the case of the open 
Hubbard and Bariev chains [35], we find — substituting (33) into (6) — that 
the eigenvectors of the transfer matrices are generated only by two classes of 
creation fields. The first class consists of the non-commuting vectors B a (X) 
satisfying the commutation relations 



B a (X) ® B b (u 
1 

Wi(X — v) 
Ww{X + v 
w 3 (X + v) 



B c {v) ® B d {X) + Wl0 ^ + V ] rjF(y)(l ® A(y)) 



■ffF(X)(I®A(v)) + 



w 3 (X + v) 
w 7 (X + v) 



r(X - v) 



w 3 (X + v) 



wio(A 



w 7 {X — v) 



F[X)B[v) 



w 6 (X - v)w 7 {X + w 10 (X - u)w 5 (X -v)~. 6 
Wl (X-v)w 7 (X-v) F{U)B{X) 



(47) 



where the functions Wi(fi) are those defined in (2) with x = e A . The second 
class of creation fields contains F(X), of (32), which commute among them- 
selves, i.e., F(X), F{y) = 0. Here, fj is a vector defined by fj = (0, 1, 1,0), I 
is a 2 x 2 identity matrix, and A denotes a submatrix given as 

\ 



A n (u) A 12 (u) 
A 2 i(v) A 22 (v) J 



(48) 



Performing the standard procedure, which is to keep the diagonal fields always 
on the right-hand sides in the commutation relations, and after several steps 
of substitutions, we arrive at the following commutation relations 



+u.t, (50) 
i / \ \ t~. / \ w\{v + A) — w 2 {X + v)wAX + v) r /, 

A^ A ^hA — ' rxrr A a. r E {r(x + u-2i 1 ) 

w^is - Xjwziv + A u>i A + v) L 



ca 

gf 



d,e,f,g=l 
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r(X-u)f b B e (u)A gd (X)}+u.t. 



(51) 



Here a,b,c = 1,2. In the commutation relations (49)-(51), we omit all un- 
wanted terms (u.t.) because they consist of a complex mixture of creation and 
annihilation fields and need a lot of space to display. The complexity of the 
unwanted terms plus the appearance of negative pseudo-momenta makes it 
very hard to perform the algebraic Bethe ansatz in a systematic way, in con- 
trast to the case of the ID Hubbard model with periodic BC [34]. However, 
we notice that the first term in each of the commutation relations (49)-(51) 
contribute to the eigenvalues of the transfer matrix which should be analytic 
functions of the spectral parameter A. Consequently, the residues at singular 
points must vanish. This yields the Bethe-ansatz equations which in turn as- 
sure the cancellation of the unwanted terms in the eigenvalues of the transfer 
matrix. Hence we prefer to use the analytical properties rather than an anal- 
ysis of the unwanted terms to derive the Bethe-ansatz equations. Fortunately, 
Eq. (51) reveals a hidden SU(l)-symmetry structure of the nesting transfer 
matrix, which is realized by the auxiliary r-matrix given by 



r(A) = 



fl 













Qi (A) 


&(A) 








f?2(A) 


f?i(A) 





vo 








1/ 



(52) 



with 



pi (A) = w 8 {\) ttt , £ 2 (A)=w 9 (A) — . (53) 



This hidden symmetry, leading to a factorization of the spin sector, plays a 
crucial role in the exact solution of the model (23). For q = e 17 , < 7 < tt, the 
matrix (52) is nothing but the scattering matrix of the anisotropic Heisenberg 
XXZ model with Boltzmann weights 

sinh2i7 sinhA 
ft(A) = sinh(A + 2i 7 )' ^ (A) = sinh(A + 2i 7 )- (54) 



Thus, the deformation parameter 7 also plays the role of an anisotropy pa- 
rameter in the hidden XXZ open chain. We also note that the commutation 
relation (47) exhibits an important symmetry, i.e., 

ft r (\\ - K(A)w 7 (A) +wi (A)w 5 (A)]w 7 (-A) ^ 
w 7 (\)w 10 (-\) 



leading to a symmetrization of the multi-particle states. 
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Following the argument of Refs. [34,35], we phenomenologically construct the 
n-particle state, 



• • • , O) = • • • , ^^'-'""lO). (56) 



Here _F a i>---> a "- are the coefficients of arbitrary linear combinations of the vec- 
tors reflecting the 'spin' degrees of freedom with Oj = 1, 2 (i — 1, . . . , n). The 
n-particle vector <3> n satisfies the recursion 



$ n (v u ...,is n ) = B ei (v l ) (g>$ n _ 1 (z/ 2 ,...,z/„) (57) 

n 

~ { [V® F( V l) • • • , Vj-i, V j+1 , ...,U n ) 

3=2 

[B^GflM, ...,!/„)- (I ® Mvj)) • ..,!/»)]} , 

where the indices = 1, 2 (i = 1, . . . , n) have been suppressed on the left-hand 
side for brevity. We remark that ff excludes the possibility of two up- or two 
down-spin particles residing at the same site. F creates a local particle pair 
with opposite spins. The coefficients Gj™\ and in turn can be determined 
from the symmetry of the wave functions 



= — - ( v $n(fl, • • • , "j+l, Vj, ■■■,Vn)- r{yj - (58) 

and the constraint arising from the cancellations of the unwanted terms in 
the corresponding eigenvalues of the transfer matrix. Explicitly, the one- and 
two-particle vectors read 



$i(z/i) = J B ei (z/ 1 ), 

$ 2 K v2) = K(yi) ® B e2 {v 2 ) + Wl °( Vl + " 2 U fM fi ® ^2) 

W3^i + V 2 ) L 

tu 7 (i/i + v 2 )w 1Q (v 1 - V 2 ) 



w 3 (^i + ^2)^7(^1 - "2) 



fjF{u{)B{v 2 ), 



(59) 



(60) 



respectively. Again, spin indices e\ and e 2 are assumed implicitly on the left- 
hand side such that (56) is fulfilled. 

According to the algebraic Bethe ansatz, the requirement that the unwanted 
terms in the eigenvalues of the transfer matrix cancel exactly yields the so- 
called Bethe-ansatz equations, which are quantization conditions for the ra- 
pidities. The property that the eigenvalues of the transfer matrix should have 
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no poles suggests an alternative way to derive the constraints on the rapidi- 
ties. It turns out that the Bethe-ansatz equations obtained by the latter way 
indeed imply the integrability of the model as can be seen by checking the 
consistency of the two-level transfer matrices. Letting the diagonal fields act 
on the state (56) and using the commutation relations (49)-(51), we get 



" sinh |(A — Vj + 2i 7 ) sinh |(A + vA . , .. . , 

= ^ 5 smhllA-^smh^A^^^ 1 ^^' " " " ' "»» + ^ ^ 

D(\)\$ n (vi, . . . ,v n )) 

^ N " cosh |(A — Vj) cosh i(A + v« — 2i 7 ) , . ,. . . 

=£ - (A) n cosh |(A-^ + 2i 7 ) cosh KA + ^ ^ +»■*■•««) 

ioo(A)|$„(vi,...,I/„)) 

" sinh \{\- Vj + 2i7) sinh \{\ + za,) 
~ £a+l( } /i sinh |(A - ^) sinh p + ^ - 2i 7 ) 

-(A + ^ - 2i 7 )^r(A - i^)Jf r(A + ^ 2 - 2i 7 )Jgr(A - 

,( A + „ B _ 2i 7 )^"V(A - ^)t£_J$n(^i, ...,!/„)> + u.t. 



xr( 
x . . .r( 



(63) 



In the last equation above, the summation convention is implied for the re- 
peated indices dj,fj and gj except for the indices a and ej. The n-particle 
state with indices tj on the right-hand side, which is as defined in (56), should 
be equivalent to the one with indices a,j on the left-hand side. It is easily 
found that (63), the eigenvalue of the submatrix A aa (A), involves a nesting 
double-row transfer matrix consisting of the inhomogeneous Lax operators 
r(A + z/ n — 2i 7 ) and r(X — u n ). For convenience, we shift the rapidities, A = w+i 7 
and Vj = Vj + i 7 , and obtain 



r(«)|$„(ui, . . .,«„)) = A(«, {vj})\$ n (vi, . . .,«„)) 

sinh \ {u — Vj + 2i 7 ) sinh \{u + Vj + 2i 7 ) 



W^ + (u + i7)ei(u + i7) [] 

j'=i 

+W 2 + (u + i 7 ) £a+1 (u + i 7 ) 



sinh |(u — f j) sinh \{u + Vj) 



x 



tt sinh \i u - + 2i 7) sinh |(m + ^ + 2i 7 ) (1) 
11 sinh Ku-Uj) sinh \{u + Vj) W>Wt) 



3=1 



+W / " 4 + (M + i 7 )cr 4 (M + i 7 ) J] 
x|$ n (ui,...,u n )) 



cosh — f j) cosh |(w + f j) 
L x cosh \{u — Vj + 2i 7 ) cosh ±(u + t> j + 2i 7 ) 

(64) 
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provided that 



Wf(u + i7)WT(u + i-i)w\ L (u + i7) 



W£(u + i7)W 2 ~(M + i-i)wl L (u + i7) 



A (1 \u = v v { Vj }), (65) 



for j = 1, . . . , n. Here A^(u, {vj}) are the eigenvalues of the nesting transfer 
matrix (67): 



r (1) (w, { Vj })F ei ' - > en = A (1) (w, { Vj })F 



ei,...,e n 



(66) 



The nesting transfer matrix reads 



r (1) (u, {vj}) = Tr c 



T (1) (m)T (1) ( 



(67) 



where 



rW( u ) 
f«( M ) 



e n gn-l 



: ri2(« + ui)};^...ri 2 (« + u„)^ Sf , 



T (1) " 



\dnfn 
V n)d, 



ian ...r 21 (u-v 1 )t{ 1 - 



(68) 
(69) 



In the previous two expressions, we used the standard notation r 12 (u) = p-r{u). 
Here p is a standard permutation operator, which can be represented by a 
2 2 x 2 2 matrix, i.e., p a /3,-yS = ^a<5<Vr We also note that r± 2 = r 2 \ for the r-matrix 
(52) and the trace operation in (67) leads to the identification g n — d n — a. It 
can be seen that the coefficients F ei '-' en act as the multi-particle vectors for 
the inhomogeneous transfer matrix (67), which characterizes the spin sector 
of the model. 



So far, we managed to solve the charge degrees of freedom. The next task, the 
diagonalization of the anisotropic Heisenberg XXZ model with open bound- 
aries, was done previously [13]. Thus we immediately obtain the eigenvalue of 
the nested transfer matrix (67) as 



A (1) (m, {itx, . . . ,u M }, {vi, ■ ■ -,v n }) 

2 sinh(w + 2^) cosh-u -q sinh(-u — Ui — 2vy) sinh(-u + u{) 



sinh 2{u + i7) fj[ sinh(-u — u{) sinh(-u + Ui + 2^) 
2 cosh(w + 2vy) sinhw sinh(w — Vj) sinh(w + Vj 



sinh 2 (u + i'y) jj[ sinh(w — Vj + 2i^y) sinh(w + Vj + 2'^) 

jl sinh(-u — ui + 2i7) sinh(w + m + Ay)) , , 

jS x sinh(-u — Ui) sinh(w + m + 2i^y) 



provided that 
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cosh (u k + 2i7) , r sinh(-u fe — i>,) sinh(w fc + Vj) 



n 



cosh 2 -u fc ^ sinh(w fe — ^ + 2i7) sinh(w fc + ^ + 2i7) 

= fi sinhK-^-2i 7 )sinhK + ^) fc = 1 M (n) 

^ sinh(« fc - + 217) sinh(u fc + w z + 417) 

l^k 

where M is the number of itinerant electrons with spin down, and n is the 
total number of itinerant electrons. The eigenvalue (70) and the constraint 
(71) on the spin rapidities u k and ui have paved the way to diagonalize the 
transfer matrix (28) completely. Making a further shift of the spin rapidities, 
i.e., Uk = Uk — 17 an d ui — fj,i — ij, the eigenvalues of the transfer matrix (64) 

r(u) |$ n (ui, . . . , v n )) = A(u, {/it}, {vj})\^ n (v u . . . , v n )) (72) 
are given by 



a/ r t r tt sirm H u ~ v i + 2i 7) siim l( M + v i + 2i 7) 

A(«, {Uzl, {f,}) = TT 2 — j-7 ^ -4t ^ - 

V 1PJ ' 1J/; sinhi(M-t; j )sinh|(M + t; i ) 

x + i7)Wf(it + i7)w 2L (w + i7) 

+W^(u + i>y)W 2 (u + \^)wl L {u + i7) 
il sinh(-u — Hi — vy) sinh(ti + fii — ij) 1 
jLi sinh(w - Hi + vy) sinh(w + /i; + i7) J 
cosh |(w — fj) cosh \ {u + Vj) 



x 



+n 



j=i cosh \ {u — Vj + 2vy) cosh |(u + + 2vy) 

x jw 4 + (u + i7)W 4 "(M + i-i)w 2 7 L {u + i7) 

+W 2 + (u + vy)W 2 ~(u + i7)w 2i (M + i7) 
^ sinh(M - /i; + 3i7) sinh(M + /i; + 3i7) j ^ 
jLi sinh(w - /ii + vy) sinh(-u + Hi + i7) J 

The rapidities of the charge and spin degrees of freedom satisfy the following 
Bethe-ansatz equations 



smh 5(^ + 17) 

-q sinhfa- - m - i7) sinhfa + m - i7) ^ 



jLi sinh(u,- - Hi + ij) smh(vj + Hi + i7) ' 
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cosh 2 (/i fc + 17) ™ sinh(/i fc - vj - 17) sinh(/i fc + ^ - 17) 
cosh 2 (/i fc - i7) sinh(/i fc - vj + 17) sinh(/z fc + vj + 17) 

vI sinh(/i fc - in - 217) sinh(/Xfc + /i/ - 217) 
sinh(/x fc - /i; + 2i7) sinh(/x fc + /!/ + 2i7) 



-q »mh(/x fc - IM, ~ 217) sinh(/x fc + ^ - 217) ^ 

1=1 
i+k 



for j = l,...,n and k = 1, . . . , M, respectively. Here, we introduced the 
notation 

C (y. e ± ) = COShl 2^-Z ± ) m 



From (22) and (73), we obtain the eigenvalues of the Hamiltonian (23) as 



" 2sin 2 7 2 

h,= } h 2L cos 7 sin 7 

j— [ cos 7 — cosh Vj cos 7 

1 , , , I, , 1 



cot-( 7 -n 



+ tan -(7 + t) + cot -(7 - r ) + tan -(7 + £") 



(77) 



Apart from a shift of the boundary parameters the Bethe-ansatz equations 
(74) and (75) coincide with those obtained by the coordinate Bethe ansatz in 
Ref. [32]. Notice that, besides the obtained Bethe-ansatz equations, we in ad- 
dition presented a systematic way to formulate the algebraic Bethe ansatz for 
Hubbard-like models with open boundary fields and obtained the eigenvalue 
of the transfer matrix (73), which is essential for the investigation of finite- 
temperature properties of the model [33] . Furthermore, the two-level transfer 
matrices, characterizing the charge and spin sectors separately, allow us to 
embed different impurities into the system. From the Bethe solutions (74), it 
is found that the boundary fields characterized by C(' u i?^ ± ) ac t indeed non- 
trivially on the densities of roots for spin rapidity fii and charge rapidity Vj, 
and thus change the ground state properties as well as the low-lying energy 
spectrum. The function C( -y i^ ± ) contributes a phase shift to the density of 
roots of the rapidities. Though the first factor on the left-hand side of the 
Bethe-ansatz equation (75) originates from the pure boundary effect of the 
spin sector, it contributes to both charge and spin rapidities in a nontrivial 
way. Of course, we may treat other boundary conditions for the model in a 
similar way. In the following section, we are going to discuss the boundary and 
impurity effects on the ground-state properties of the model at half filling. 
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4 The ground-state properties 

Due to the boundary effects, the Bethe-ansatz equations (74) and (75) are 
not merely a doubling of the Bethe equations for the periodic Ujosp(2|2)] 
chain [26]. Thus the boundary fields contribute nontrivially to the ground- 
state properties of the model. If we consider the ground state at half filling, 
corresponding to the case that the Vj are real, while the spin variables form 
strings of type /i + i|, the discrete Bethe-ansatz equations (74) and (75) may 
be written as 

2L9 1 (^v v ± 7 ) - 6 2 (\vj, |n - 6 2 (\vj, ±r) 

M 

= 2nlj - ]T [9 2 ( Vj - fi h 7 ) + 9 2 ( Vj + fi h 7)] , (78) 
i=i 

2^1(^,7) 

M 

= 2nJ k - [Oifak ~ Hi, 27) + Oi(n k + in, 27)] 

l=i 
l+k 

n 

~ S [^(^ ~ v ii 7) + 2{fik + Vj, 7)] , (79) 
i=i 

for j = 1, . . . , n, and k — 1, . . . , M, respectively. Here we introduced the shift 
functions 

2 arctan (tanh v cot 7) , (80) 

2 arctan (tanhf tan 7) . (81) 

The integers Ij and may be regarded as quantum numbers associated to 
the Bethe-ansatz equations. If we define I-j = —Ij, J^t = —Jk, v-j = —Vj, 
and ji-i = —fii, and pass to the thermodynamic limit L — > 00, n — > 00, and 
M — > 00, with n/L and M/L kept finite, the thermodynamic Bethe-ansatz 
equations read 

+ 2 / d/i— ^(t;-/i,7)p s 00 (^) , (82) 

1 f 1 d °r d 

P S co(^) = 2^ ^d^^(/x,7)+ y du— 2 (Ai-u,7)p»(u) 



0i (u, 7) = 

^2 (U, 7) = 



— i In 



-i In 



sinh(i7 



sinh(i7 + v) 
cosh(i7 + v) 
cosh(i7 — v) 
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+ / dv— 1 (ii-v,2i) / f oo (v) . (83) 

oo J 

Here p^ (f) and p^p) denote the densities of roots of the charge and spin 
rapidities at half filling, respectively. The term 

M«, t±) = -02(v, 7 ) - e 2 a v , k+) - e 2 a v , §r ) (84) 



characterizes the charge contributions to the densities of roots from the bound- 
ary potentials, whereas 

sb (/x, 7) = 20x (p, 7 ) - 9 1 Gu, 2 7 ) - 0x (2/i, 2 7 ) - 2 (/x, 7) (85) 



denotes the boundary contributions in the spin sector. We would like to stress 
that, although c b(f , ^ ± ) arises completely from the charge degrees of freedom, 
and s b(p? 7 ) only from the spin degrees of freedom, both terms contribute 
nontrivially to the densities of roots p^ (f ) and p^,(/x). We also find that the 
ground state of the system is a singlet, which means that the variables ji and 
v occupy the entire interval from — oo to oo. By Fourier transformation, the 
solutions of (82) and (83) have the form 



ploiv) = ^ / ^ v d ^ = pSH + j/bM, (86) 

— oo 

pLM = ^ / £oM e_i ^ = PoM + 7/bM- (87) 

— oo 

Here, Po(o;) an d PoO^) denote the bulk densities of roots for the charge and 
the spin rapidities, which are given by 

2sinh(f - 7 )o; 1 

PoM " coshfc ' PoM - c^hf? (88) 



respectively. The remaining parts p\{oj) and p\{uS) contain the contributions 
caused by the boundary terms, i.e., 



PbM = 
pIH = 



[sinh £ + cj + sinh £ cj] cosh(| — 7 )u; sinh ^a; 



sinh(7r — 7 )cj cosh |a; 
[sinh t; + u + sinh £ _ u;] 



sinh(f - 



+ 



sinh(| — 7 )u; 



2 sinh(7r — 7 )cj cosh |a; 2 sinh(| 



cosh 



,-u; 



(89) 
(90) 
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respectively. In the expressions (89) and (90), we separated the effects of the 
boundary potentials and the pure boundary effects, i.e., corresponding to the 
first and the second terms on the right-hand sides of (89) and (90), respectively. 
It is straightforward to recover the result for free BC for the model (23) by 
switching off the boundary potentials via ^ — > 0. Because the open BC do 
not spoil the symmetries p%o(co) = p^ (— cj) and (^(u) = p^ c (— cu), the energy 
per site of the singlet ground state (n — -|), calculated from (77), reduces to 
the form 



Tcoshd — 7)0; sinh(7r — 7)0; 
E = -4 sin 7 / - 



cosh |u; sinh ixoj 




du + 2 cos 7 



+ sin 7 



cot-(7-£ H 



+ tan ^(7 + f +) + cot ^(7 - D + tan ^(7 + C) 



(91) 



The first term in the ground-state energy (91) is the bulk ground-state energy 
which coincides with the result of the periodic chain [26]. We emphasize that 
the boundary potentials do not only enter in the expression for the ground 
state energy explicitly as cot |(7 — ^ ± ) + tan ^(7 + £ ± ), but also implicitly via 
pUv) and pUv) of (86) and (87). 

Before closing this section, we discuss the problem of embedding integrable 
impurities into the open chain [20-22] . The algebraic Bethe ansatz provides us 
with a natural way to incorporate different kinds of impurities. If we embed 
two impurity vertices at the boundaries, namely, extend (20) to 

T(A) = R rfi (X + Pt )Rl,o( x ) r l-i,o( x ) ■ ■ ■ R2fl(X)R lfi (X)Rt fi (X + p t ), (92) 



the impurity Hamiltonian with the impurity-host charge interactions and ex- 
change coupling between the impurities and boundaries can be determined 
by 



{Str [^ + (0) J R; i0 ( Pr ) J R r " 1 ( Pr )] 



Str K + (0) 

+Str [^ + (0) J R r ,o(p r ) J R , L , (0) J RZ, 1 o(°) i? r;o 1 (Pr)]} 

+ ^R lfi (0)R e , ( Pe )KU0)Rii( P i)R^(0) 
+R w (0)R' etO ( Pe )R^(p e )R^(0). 



(93) 
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112 L-1 L r 

Fig. 1. Impurities coupled to each of the boundaries. 

Due to rather lengthy algebra, we schematically present the interactions be- 
tween boundaries and impurities by Figure 1 instead of presenting ex- 
plicitly in terms of fermionic operators. The quantities p T and pe shifting the 
pseudo-momenta characterize the impurity rapidities or impurity strength. As 
discussed previously [11,21], this shift preserves the integrability of the model 
and allows one to continuously vary the strength of the impurity coupling 
to the boundaries revealing different impurity effects. This embedding of the 
impurities, carrying both charge and spin degrees of freedom, leads to an ad- 
ditional factor 

tt sinh \{vj +p m - 17) sinh \{vj - Pm - 17) , Q4 \ 
J-Jl r sinh \ (vj +p m + vy) sinh \ (vj -pm + i"/)' 

on the left-hand side of the Bethe equation (74), thus contributing to the 
densities of roots for the charge and the spin rapidities in (86) and (87) at 
order 1/L, namely, 



pU^=Po^) + j[pI^)+pK^], (95) 
=ft(u) + j\ft(u)+ft(u)]. (96) 
The contributions to the densities of roots from the impurity terms are 



p\{uj) = (cosp r u + cos p e u) p c (uj), (97) 
pKu) = (cos p r uj + cos peuj) p s (uj). (98) 

In this case, the ground state energy (91) has to be changed slightly by 
[pl(uj) + (to)} replacing pX(u), which affects the surface energy and finite- 
size corrections. 

From the above discussion, one can easily distinguish the effects of the im- 
purities, boundary potentials, and the free edges of the system, i.e., the ef- 
fects of dynamic magnetic impurities, the external scalar boundary fields, and 
open boundary conditions. In addition, the integrability of the open chain 
also allows us to embed a forward scattering impurity (without any reflection 
scattering amplitude) into the bulk part. Although the resulting impurity 
Hamiltonian is different from the boundary Hamiltonian, the Bethe ansatz 
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shows that the impurity effects do not depend on the position of the impu- 
rity due to the pure forward scattering that is required by integrability [11]. 
Another interesting embedding of the impurity would be given by operator- 
valued boundary i^-matrices [22] which lead to Kondo-impurity-like terms 
in the Hamiltonian. These different embeddings of the impurity would provide 
results essential to the study of the thermodynamic properties such as low- 
lying excitations, finite-size corrections, magnetization, etc. The other string 
solutions to the Bethe equations (74) and (75), which form the charge and 
spin bound states, give an independent approach to the investigation of the 
low-lying excitations for the model. We intend to consider this situation in the 
future. 



5 Conclusion and discussion 



In summary, we have discussed the algebraic Bethe-ansatz solution for the 
extended Hubbard model with boundary fields associated with the quantum 
superalgebra U g [osp(2|2)] in terms of the graded QISM. Two classes of solu- 
tions of the graded RE leading to four kinds of possible boundary terms in the 
Hamiltonian were obtained. The Bethe-ansatz equations, the eigenvalue of the 
transfer matrix and the energy spectrum were given explicitly. The ground- 
state properties in the thermodynamic limit were also studied. We found that 
the model exhibits an anisotropic Heisenberg XXZ open chain as its nesting 
transfer matrix characterizing the spin sector. This nesting structure seems 
to be different from that of other extended Hubbard models [22] . The Bethe- 
ansatz results allow us to embed impurities at the boundaries of the model. 
Our results provide a useful starting point for studying the thermodynamic 
properties and correlation functions for the model. The boundary potentials, 
pure boundary effects and impurity effects contribute nontrivially and sepa- 
rately to both the density of roots of the charge rapidity and the spin rapidity 
at order 1/L. The impurity strengths pi and p r , the boundary parameters ^ 
and the g-deformation parameter 7 change the asymptotic behavior of the 
thermodynamic Bethe-ansatz equations, i.e., the band filling, magnetization, 
susceptibility, compressibility, finite-size corrections, etc. Our computations 
could open an alternative way to study the thermodynamic properties for the 
integrable double sine-Gordon model [30] as well as to investigate the tunnel- 
ing effects in quantum wires [31]. 
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A The boundary .R^-matrices 



Let us consider an algebraic ansatz which fixes the boundary /^--matrices. 
Due to the isomorphism between K + - and fT_-matrices, we need to solve 
the RE (6) to determine them. Let us first fix the left boundary i^„-matrix 
(11). By substituting it into the RE (6), one may find that the RE (6) in- 
volves two variables A and v which make the functional equations involving 
K { 1\\), i = 1,...,4 much more complicated. Nevertheless, taking into ac- 
count the structure of the i?-matrix (1) and the RE (6), we find that the 
functional equations of the RE arising from the positions corresponding to 
the permutation operator P a p,-yS — ( — l) P ^ P ^^a<s^37 provide us with simpler 
equations than the ones arising from other positions in the RE. These simpler 
equations allow us to separate K^(X), K^{y) into factorized forms as follows 





[v) 


w 3 (X + u)w 2 (X - 




(A) + w 3 (X - 


v)w 4 (X + 




(A) 


K m 


V) 


w 2 (X + v)w 3 (\ - 




(A) +w 4 (X - 


v)w 3 (X + 


v)K ( ^ 


(A) 


K w 


[v) 


w 3 (X + u)w 2 (X - 




(A) + w 3 (X - 


is)w 4 (X + 


u)K (3) 


(A) 


K m 


V) 


w 2 (X + v)w 3 (X - 


v)K^ 


(A) +w 4 (A - 


v)w 3 (X + 


u)K^ 


(A) 


K w 


V) 


w 3 (X + u)w 2 (X - 


v)K® 


(A) + w 3 (X - 


v)w 4 (X + 




(A) 


K w 


V) 


w 2 (\ + u)w 3 (X - 




(A) + u> 4 (A - 


v)w 3 (X + 




(A) 



kk.2) 



U.3) 



Substituting the Boltzmann weights of the i? 12 -matrix (1) into the equations 
above, and analyzing the structure of these equations, we can infer that the 
functions on the right-hand side of each of the equations (A.1)-(A.3) involving 
the variable A should cancel because the left-hand sides of (A.1)-(A.3) are 
functions of the variable v only. This cancellation leaves us with the following 
expressions for K^(X) 



K { y (A) = {x + C^) (x + C^) (x + C^) , (A.4) 
K {2 \X) = (x- 1 + CW) (x + C™) (x + C®) , (A.5) 
xi 3) (A) = {x + C™) (x- 1 + C™) (x + C®) , (A.6) 
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K^(X) = (x + CW) (V 1 + C^) (x~ l + C®) , 



(A.7) 



with a minimal number of coefficients that are to be determined. We regard 
Eq. (A.4)-(A.7) as the main structure of the i^-matrix. Here x = e A as in 
section 2. Employing the RE again with (A.4)-(A.7) it is easily found that 
only one coefficient is free, and we have the relations = = —C^q 2 
or C« = -C ( Vq 2 = -q 2 /C (2 \ The first case yields the solution (13) while 
the second case corresponds to (14). The corresponding boundary terms have 
been given by (25)-(27), respectively. Due to the fact that the left and right 
boundaries are independent of each other, the two cases of boundary terms 
allow four possible combinations. In particular, we find a spin-degenerate sit- 
uation with chemical potentials at the boundaries (i) identical for spin-up and 
spin-down electrons, of (26), or (ii) different yielding (27). Or, the spins are 
distinguished either (iii) at the right boundary, i.e., Vl,± as i n (27), or (iv) 
at the left boundary, i.e., V^± as in (27), and other chemical potentials given 
by (26). In all cases, U(£±) is given by (25). These potentials constitute the 
general integrable boundary terms corresponding to the diagonal boundary 
X-matrices. In this paper, we restrict our discussion to the physically most 
realizable, symmetric case given by (13). 



B The Hamiltonian 



Because of Str [-ft'+(0)] = 0, we have to consider the second-order expansion 
of the transfer matrix (19) with respect to the spectral parameter A in order 
to construct the Hamiltonian with boundary fields. Following [15], we derive 



L-l 



H= £ H jJ+1 + -K'_(0) + { Str \K' + (0)R' LO (0)P LO 

3=1 ^ 



+ 2Str 



K + (0)R" LO (0)P LO ] + -Str [k + (0)(R' L0 (0)P L0 ) 2 \/ 



Str <(0) +2Str K + (0)R' LO (0)P LO 



(B.l) 



where 



H 



dR 



dA 



(A) 



(B.2) 



A=0 



The prime denotes the derivative with respect to the spectral parameter A. 
Using (B.2), we can write the bulk Hamiltonian in terms of two commuting 
species of Pauli matrices a and r, i.e., 
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V|r/ + aJ +1 r/ +1 l + r/r/ +1 ' 



r i r i+i + r j r i+i 



4 2{q-q- 1 )) 
4 2(g-(z-i). 



g + g 1 

9 + 9 _1 
~8(9-9- 1 ) 

+l(^- (7 l + i + r /- 7 i + i) 



a*r* +1 + ^ +1 r/ + o]r] + ^ +1 r/ +1 + 4 1 



(B.3) 



For the solution (13) of the RE, the second term in (B.l) gives the boundary 
terms at site j — 1 



;i + g^-alT? - (29 - + 9 2 £-)(^ + rf 
4(g-£_)(l + g£_) 



(B-4) 



The boundary terms at site j — L are then given by the remaining terms of 
(B.l) as 

^ = 4(9-e + )(l + ^ + ) • (R5) 



Using the Jordan- Wigner transformation [36], apart from a factor —1/(9— q~ v ) 
- absorbed in the constant c 2 in (22) — the bulk Hamiltonian (B.3) together 
with the boundary terms (B.4) and (B.5) has the form presented in Eqs. (23)- 
(26). The last term in (B.3) should be taken into account in the boundary 
terms. 



C Useful commutation relations 



For the factorization of the transfer matrix (19) acting on the pseudo- vacuum 
state, we need the following commutation relations 



w l {2X)C 1 {X)B 1 {X)=w 2 {2X) [B(X)B(X) - A 11 (X)A 11 (X) 
w 1 {2X)C 1 (X)B 1 {X) =w 4 {2X) [B(X)B(X) - A 11 (X)A 11 (X) 
Wl (2X)C 2 (X)B 2 (X)=w 2 (2X) [B(X)B(X) - A 22 (X)A 22 (X) 
w 1 (2X)C 2 (X)B 2 (X)=w 4 (2X) \B(X)B(X) - A 22 (X)A 22 (X) 



(C.l) 
(C.2) 
(C.3) 
(C.4) 
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w 1 (2X)C 3 (X)F(X) = w 5 (2X) [B(X)B(X) - D(X)D(X) 

-w 2 (2X) [C 4 (A)^(A) + C 5 (X)E 2 (X)] , (C.5) 
w 4 (2X)C 3 (X)F(X) = w 2 (2X) [A U (X)A 11 (X) - D{X)D(X)] - C 4 (X)E 1 (X) 

+w 5 (2X)C 1 (X)B 1 (X) - w 8 (2X)C 5 (X)E 2 (X), (C.6) 
w 4 (2X)C 3 (X)F(X) =w 2 (2X) [A 22 (X)A 22 (X) - D(X)D(X)] - C 5 (X)E 2 (X) 

+w 5 (2X)C 2 (X)B 2 (X) - w 8 (2X)C 4 (X)E 1 (X), (C.7) 

which can be derived directly from (34). From these relations, we obtain 



C 3 (A)F(A) 



(? 2 " 1) 



(x 2 q 2 — l)(x 2 — q) 



2(q 2 - 1) 



(x 2 + l)(x 2 - q 2 ) 



B(X)B(X) + 



A 11 (X)A 11 (X) + A 22 (X)A 22 (X) 
2(q 2 - 1) 



{x 2 q 2 - \){x 2 + 1) 



D(X)D(X), (C.8) 



C,(m(X) = -gg ]l A u (A) A u (A) 



(rrV-r)^ 2 -? 2 )' 

(xY-i)(x 2 -? 2 ) 22( } 22( } " ^V^T D(A)D(A) ' 



C 5 (A)£ 2 (A) 



x 2 (q 2 -l) 2 (q 2 + l) , 



(xV-l)(x 2 -g 2 ) 



A n (A)A n (A) 



+ ^^4i 22 (A)i 22 (A) - -4^(A)Z)(A). 



(iY-i)(s 2 -g 2 )' 



X 2g2 _ ]_ 



(C.9) 



(CIO) 



In addition, one also can show that Ci(X)Bj(X) = for i ^ j, i = 1, 2, 3, and 
j = 1,2. With the help of (C.8)-(C10), it is not difficult to derive (35). 
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